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Enhancing the signal-to-noise ratio of solar-like targets 
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The analysis of the first solar- like targets done by CoRoT has shown that the oscillation amplitudes are about 25% below 
the theoretical amplitudes while the convective backgrounds are up to three times higher than in the solar case (Michel et 
al. 2008). In such conditions, the comb-like structure of the acoustic modes has smaller signal-to-noise ratios than initially 
expected complicating the characterization of individual modes. In the present work we apply the curvelet filtering to the 
solar-like targets already observed by CoRoT as well as a partial reconstruction of the signal from the obtained spacing of 
the comb-like structure of the acoustic modes. It enables us to enhance the signal-to-noise ratio of the ridges in the Echelle 
diagrams. Finally, we study how the analysis of the p modes can be improved. 
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1 Introduction 

The signal-to-noise ratio (SNR) of solar-like stars can be 
enhanced by taking advantage of the asymptotic properties 
of their acoustic (p) modes. In the Fourier spectrum of these 
stars, the peaks that correspond to the acoustic modes are 
almost equally spaced in a given frequency range. Here we 
describe two methods: the partial reconstruction of the sig- 
nal and the Curvelet filtering and we apply them using as- 
teroseismic observations obtained by CoRoT (Convection, 
Rotation and planetary Transits). 



2 Partial reconstruction of the signal 

To increase the SNR of the periodic structures (the modes) 
we perform a selective filtering of the power spectrum of the 
power spectrum (PSPS) (see for details: Regulo and Roca 
Cortes 2002; Mathur et al. 2010a). 

We start by computing the PSPS of a short range in fre- 
quency (< 1000/iHz) around the maximum power of the p- 
mode band (the asymptotic part). The resultant PSPS is fil- 
tered by multiplying it by a window function, which is equal 
to 1 for all the equally spaced bins around multiples of the 
large separation of the p modes (that should have been esti- 
mated before) starting at zero, while the rest of the bins are 
settled to 0. The inverse Fourier Transform of this filtered 
PSPS produces a "recovered" power spectrum of the stellar 
p modes with a higher SNR compared to the original one. 



We have applied this technique to two high SNR CoRoT 
targets: HD181420 (Barban et al. 2009) and HD49933 (Ap- 
pourchaux et al. 2008; Benomar et al. 2009). An example of 
the results for the second star is shown in Fig.[T] 
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Fig. 1 Raw (top) and recovered (bottom) power spectrum 
density (PSD) of HD49933. The comb-like structure of the 
modes is enhanced in the recovered PSD. 



Fig. [21 shows the raw and the recovered PSD of two 
other low SNR CoRoT targets: HD181906 (Garcia et al. 
2009) and HD175726 (Mosser et al. 2009). In both cases, 
the modes are very difficult to be distinguished in the raw 
PSD but in the recovered ones, we can easily see the comb- 
like structure of the p modes. 
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Fig. 2 Raw (top) and recovered (bottom) PSD of the low SNR targets HD181906 and HD175726. 



This method tends to rigidify the modes (constant spac- 
ing) towards the borders of the considered region where the 
modes could be less asymptotic and the large separations 
could change significantly with frequency. The recovered 
amplitudes and line-widths of the modes are unfortunately 
modified, thus this method is only reliable, for the moment, 
to study the frequency distribution of p modes. 

3 Curvelet filtering 

It is commonly accepted that Echelle Diagrams can help in 
the process of mode tagging as it has already been done in 
early stages of helioseismology (e.g. Grec et al. 1981). To 
improve the SNR, Bedding et al. (2004) proposed to use 
this Echelle diagram by smoothing it in the vertical direc- 
tion. However, this method only works well in the frequency 
region in which the ridglets are quasi- vertical, i.e., in the 
asymptotical regime. Even more, a very good a priori of the 
large separation is needed. 

We propose here another denoising technique based on 
the use of Curvelets. Curvelet transforms are built to deal 
with curved structures of finite size in the image. A curvelet 
transform (Candes & Donoho 1999) is a ridgelet transform 
(Candes 1998) but used in a a localized manner: at a small 
scale any curved line can be approached by a straight line. 
Using standard wavelets, many coefficients are needed to 
well define the curved figures while with the curvelets, a 
few coefficients are enough. 

The curvelet transform has been successfully tested in 
asteroseismology by an extensive use of Monte-Carlo sim- 
ulations and real solar data (Lambert et al. 2006). The fil- 
tering becomes more efficient when the SNR is small and it 
only needs a rough estimation of the large separation. 

Fig. [5] shows the raw and the curvelet filtered PSD of 
HD181906 (Garcia et al. 2009). The reduction in the noise 
level of the filtered spectrum allows to retrieve the comb- 
like structure of the acoustic modes centered at 1 800 /iHz, 
with the p-mode hump visible from around 1000 to around 
2600/iHz. 
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Fig. 3 Raw (top) and curvelet filtered (bottom) PSD of the 
low SNR target HD181906. 



The Curvelet filtering also works well for short runs, as 
in the case of HD 175726 in which less than 30 days of ob- 
servations were available (see Mosser et al. 2009). In Fig.|4l 
we have plotted the original Echelle diagram and the col- 
lapsed spectrum (sum of each horizontal line of the Echelle 
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Fig. 4 Echelle diagram and collapsed spectrum computed as indicated in the text for the CoRoT short-run target 
HD 175726: raw data (left panel) and curvelet filtered data (right panel) 



diagram) for both, the raw data and the Curvelet filtered one. 
In the Curvelet filtered diagrams we can clearly identify the 
presence of two ridges corresponding to the odd and even 
modes. 

A comb-like pattern is clearly visible in the Curvelet fil- 
tered spectrum centred at 2120/iHz (see Fig. [5]). We have 
compared the main peaks in the p-mode hump with those 
obtained by the partial reconstruction method described in 
section 2. The frequencies of the peaks obtained with both 
methods coincide very well, except for the two modes at 
low frequency for which the partial reconstruction method 
maintain a regular distance while in the Curvelet filtering 
spectrum the peaks are shifted from the regular spacing. 

4 Conclusions 

In this work we have shown two powerful methods to en- 
hance the comb-like structure produced by acoustic modes 
propagating in the interior of stars. The partial reconstruc- 
tion method allows to obtain the asymptotic large spacing 
by computing the PSPS and the reconstruction of the pattern 
in the PSD. The curvelet filtering technique uses the Echelle 
diagram to filter the noise and enhance the curved ridges of 
the odd and even modes. In both cases, the identification or 



tagging of the p-mode ridges is simplified and it can poten- 
tially help in the determination of the rotation speed and the 
inclination angle of the star, parameters that are correlated 
(Gizon & Solanki 2003; Ballot, Garcia & Lambert 2006) 
and are very difficult to disentangle even with modes ob- 
served with high SNR (e.g. Appourchaux et al. 2009; De- 
heuvels et al. 2010). 

Both methods have been intensively tested using Aster- 
oFLAG simulated data (Chaplin et al. 2008) as well as on 
public CoRoT data, as shown here, and in the new low SNR 
target HD170987 (Mathur et al. 2010b) giving very promis- 
ing results. Both techniques will be incorporated to our au- 
tomatic asteroseismic pipeline (Mathur et al. 2010a) in or- 
der to analyze very low SNR targets of the survey phase of 
the Kepler mission (Bedding et al. 2010; Stello et al. 2010; 
Chaplin et al. 2010; Grigahcene et al. 2010). By doing so, 
we hope to extract a few individual frequencies in order to 
improve the determination of the stellar parameters in these 
low SNR targets (e.g. Stello et al. 2009). 
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Fig. 5 Curvelet filtered PSD of HD175726. The dotted lines are the frequencies found by Mosser et al. (2009). Blue 
dotted lines are the frequencies ui in table 4 of Mosser et al. (2009), while the red ones are those labeled as 1^2. Dashed 
green lines are the frequencies obtained by the partial reconstruction method described in section 2. 
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